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We study the dipion transitions T(4S) > YT(nS)r^-- (n = 1,2). In particular, we 
consider the effects of the two intermediate bottomoniumlike exotic states Z,(10610) and 
Z,(10650) as well as bottom meson loops. The strong pion-pion final-state interactions, 
especially including channel coupling to KK in the S-wave, are taken into account model- 
independently by using dispersion theory. Based on a nonrelativistic effective field theory 
we find that the contribution from the bottom meson loops is comparable to those from 
the chiral contact terms and the Z;-exchange terms. For the Y(4S) — Y(2S)atz~ decay, 
the result shows that including the effects of the Z,-exchange and the bottom meson loops 
can naturally reproduce the two-hump behavior of the mm mass spectra. Future angular 
distribution data are decisive for the identification of different production mechanisms. For 
the T(4S) + T(1S)r*«- decay, we show that there is a narrow dip around 1 GeV in the 
77 invariant mass distribution, caused by the final-state interactions. The distribution is 
clearly different from that in similar transitions from lower Y states, and needs to be verified 
by future data with high statistics. Also we predict the decay width and the dikaon mass 
distribution of the T(4S) > Y(1S)K* K- process. 
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I. INTRODUCTION 


The processes of dipion emission of the bottomonia T(mS) > Y(nS)z« are important for un- 
derstanding the heavy-quarkonium dynamics and low-energy QCD. Because the bottomonia are 
expected to be nonrelativistic and compact, the method of the QCD multipole expansion [1—4] 
is often used to study these transitions, where the pions emitted come from the hadronization of 
soft gluons. Though successful in describing many Y(mS) — Y(nS)rz« processes, a well-known 
anomaly about the method of the QCD multipole expansion is that it cannot reproduce the two- 
hump behavior in the experimental mm invariant mass spectra of the decays T(39) — Y(1S)a7 
and T(4S) > Y(2S)r*«- [5-7]. In a previous study [8], we found that by including the effects of 
the two bottomoniumlike exotic states Z,(10610) and Z,(10650) discovered by the Belle Collabo- 
ration [9, 10] as well as the 77 final-state interaction (FSI), the anomaly of the T(3S) — Y(1S)a7 
process can be naturally explained. Such an analysis is a modern version of the much earlier studies 
in Refs. |11, 12], where an isovector bbqğ state was considered. Although the direct decay of Z, into 
'Y (4S) is kinematically impossible, it may be illuminating to analyze the effect of the Z,-exchange 
mechanism in the T(4S) > Y(1S,2S)x^7 processes, which is performed in this study. In this 
context it is important to note that improved data on Y (nS) decays are to be expected from Belle- 
II that will start operation soon—for a detailed discussion of prospects for various decays relevant 
for this study we refer to Ref. [13]. 

The YT(4S) meson is above the BB threshold and decays predominantly to B B, so loop effects 
with intermediate bottom mesons may play an important role in T(4S) — Y(nS$)z« (n = 1,2). 
Also, the inclusion of the loops will introduce non-analyticities arising from the BB threshold 
needed to be taken into account in dispersion theory, which will be discussed later. Because 
the bottomonia are close to the open-bottom meson production threshold, the velocity of the 
intermediate bottom mesons is small and can be treated as an expansion parameter to build power- 
counting rules in a nonrelativistic effective field theory (NREFT) [14-16]. Within the NREFT 
scheme, we will calculate the dominant box diagrams in the dipion emissions of T(4S), and find 
that their contribution is comparable in size to the chiral contact terms and the Z-exchange graphs. 

In the Y (4S) — T(1S)z process, the dipion invariant mass reaches above the KK threshold, 
so the coupled-channel FSI in the S-wave is strong and needs to be taken into account. Based on 
analyticity and unitarity, dispersion theory can achieve this in a model-independent way. In this 
study, we will use dispersion theory in the form of modified Omnès solutions, in which the left- 


hand-cut contribution is approximated by the sum of the Z,-exchange mechanism and the bottom 
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meson loops. At low energies, the amplitude should agree with the leading chiral results, so the 
subtraction functions can be determined by matching to chiral contact terms. For the leading 
contact couplings of two S-wave bottomonia to an even number of light pseudoscalar mesons, 
we will adopt the Lagrangian given in Ref. [17], constructed in the spirit of the chiral and the 
heavy-quark nonrelativistic expansions. 

This paper is organized as follows. In Sec. II, we present the theoretical framework and elaborate 
on the calculation of the amplitudes as well as the dispersive treatment of the FSI. In Sec. III, we 
fit the experimental data of the v7 invariant mass distribution to determine the coupling constants, 


and discuss the contributions of different mechanisms. A summary will be given in Sec. IV. 


II. THEORETICAL FRAMEWORK 


A. Lagrangians 


Because in the heavy-quark limit the spin of the heavy quarks decouples, it is convenient to in- 
troduce the heavy quarkonia and heavy hadrons in terms of spin multiplets. One has J = Y -o +r, 
where Y and 7 annihilate the Y and 7 states, respectively, and ø contains the Pauli matrices [18]. 
The bottom mesons are collected in Ha = Va: o + P, with P,(Va) = (B(9-, BE), BO) and 
H, =—-V,-o+ P, with P,(V,) = (B+, B®, BY) [19]. 

The effective Lagrangian for the contact Y'Y^rz and YY’KK coupling, at the lowest order in 


the chiral as well as the heavy-quark expansion, reads [8, 17] 
-jiy py iy ay 
Mao = pud J')(uyu") + v J')(uyu,)v"v" + h.c. , (1) 


where v = (1,0) is the velocity of the heavy quark. The Goldstone bosons of the spontaneous 


breaking of chiral symmetry can be parametrized as 


Up, — i (utdvu — udu) ; u = exp (+s) ; 
ar + Zi mt Kt 
b=] Attn KO |. (2) 
K- K? -hs 


where F is the pseudo-Goldstone boson decay constant, and we will use Fp = 92.2 MeV for the 
pions and Fg = 113.0 MeV for the kaons. The two operators in Eq. (1) both scale as O(q2) in the 


expansion in (soft) pion momenta q;. 
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The leading Lagrangian for the Z,Y' interaction, which is needed in the calculation of the 
mechanism Y(mS) > Zyn — Y (nS)nn, reads [18]! 
Larn = X, X Catus Y (18)(Zi; upju” + hic. , (3) 
j=1,2 1 
where Zp; and Zz are used to refer to Z,(10610) and Z(10650), respectively. The Zp states are 


collected in the matrix as 


d Zi 0 


v 
Z=] ZA -52% o0]. (4) 
0 0 0 


Note that since strange partners of the Z states, Zs, have not been observed, we set the corre- 
sponding matrix entries in Eq. (4) to zero. 
'To calculate the box diagrams, we need the Lagrangian for the coupling of the S-wave bottomo- 


nium fields to the bottom and antibottom mesons [14], 
"T 
Lynn = SEF (tHe JOH.) + hc., (5) 


= 
where AOB = A(Ö B) — (d A)B. We also need the Lagrangian for the axial coupling of the 
Goldstone bosons to the bottom and antibottom mesons, which at leading order in heavy-flavor 


chiral perturbation theory is given by [19-28] 
— In at Bi iiri 
Jugo = z Hao . Uab Hp) = punti $ Uba), (6) 
where u! = —\/20'6/F + O(4?) corresponds to the three-vector components of u, as defined in 


Eq. (2). Here we will use gr = 0.5 from a recent lattice QCD calculation [24].? 


B. Power counting of the loops 


Since the T (4S) is above the BB threshold and decays predominantly into BB pairs, the loop 
mechanism with intermediate bottom mesons may play a significant role in the bottomonium 
transitions T(4S) — Y(nS)n^n-. In this section, we will analyze the power counting of different 


kinds of loops, based on NREFT [14-16]. In NREFT, the expansion parameter is the typical 


! Here we only include the terms relevant to the Y coupling rather than the full spin multiplet defined before as 


J = Y -o +m. In this way, we avoid the discussion of the internal spin structure of the Z, states, which depends 
on specific models for Z and is not really settled yet. 
? The precise value quoted in Ref. [24] is gx = 0.492 + 0.029. 
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FIG. 1: Feynman diagrams considered for the T (mS) — Y(nS)nrr processes. The crossed diagrams of (b1), 
(c1), (b2), and (c2) are not shown explicitly. The gray blob denotes the final-state interaction. 
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FIG. 2: Subleading contributions for T(mS) — Y(nS)rz that are suppressed in comparison to the four- 
point functions in Fig. 1 and hence not considered in the calculations. The corresponding power counting 


arguments are given in the main text. 


velocity of the intermediate heavy meson, namely v = 4/|myus) — Mg) — Mge |/Mge, and 


v « 1 since T(lS) are close to the B (9) B®) thresholds. Each nonrelativistic propagator is counted 


5. More details of the power counting rules are 


as 1/v?, and the full integral measure E] d^l as v 
elaborated in Ref. [15]. 

Without considering the FSI, there are five different kinds of loop contributions, namely the 
box diagrams displayed in Fig. 1 (c1), (d1), and the triangle diagrams displayed in Fig. 2 (a)-(c). 


We analyze them one by one as follows: 


1. Box diagrams, namely Fig. l(cl), (dl): As indicated in Eq. (6), the vertex for the axial 
coupling of the pion to the bottom mesons is proportional to the external momentum of the 
pion q. Both the vertices for the initial and final bottomonia are in a P-wave, and the 


product of the two vertices can be counted as O(v?), so the box diagrams are counted as 
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vq [y* = q2 [v. Note that these contributions thus have the same scaling in pion momenta 
as the leading Y Y'z contact terms from the Lagrangian Eq. (1), but are formally enhanced 


by 1/v in the non-relativistic velocity parameter. 


. Fig. 2(a): The leading B®) BC) am vertex comes from the covariant chiral derivative term 


UTI (Do us Hy) = (Hi (ig —iVo)ba Hy) [25, 26], in which the pion pair produced by the vector 
current, V^ — 3(uldHu + uO" ui), cannot form a positive-parity and C-parity state, so this 
leading vertex does not contribute to the Y (mS) — Y (nS)z7 processes. Isoscalar, PC = ++ 
pion pairs only enter in the next order O(q2) from point vertices. For the vertices of the 
initial and final bottomonia, both of them are in P-waves, so the product of them can be 
counted as O(v?). These diagrams hence count as v°v2q?/v® = vq2, and are suppressed 


compared to the contact terms « c1,» by the factor v. 


. Fig. 2 (b), (c): The leading Y (L9) B? Br vertex given by U Hi Hio, [27] is proportional 


a 


to the energy of the pion, E; ~ qr. In Fig. 2 (b), the vertex for the initial bottomonium is in 
an S-wave, and the vertex for the final bottomonium is in a P-wave, so the loop momentum 
must contract with the external momentum qr and hence the P-wave vertex scales as O(q;). 
For this reason, Fig. 2(b) is counted as v°q3/(v°mpg) = q3/(vmg), where the factor mg 
has been introduced to match the dimension with the scaling for cases 1 and 2. Analogous 
arguments hold for Fig. 2(c). This class of diagrams is therefore suppressed in the chiral 


expansion in pion momenta, compared to the c1,2 terms. 


We find thus that according to the power counting the box diagrams are dominant among the 
loop contributions, and the only ones not expected to be suppressed relative to the tree-level 
contact terms. We will therefore only calculate these in the present study. Note that all (box and 
triangle) loop contributions discussed here are ultraviolet-finite, and do not require the additional 


introduction of counterterms. 


C. Tree-level amplitudes and box diagram calculation 


First we define the Mandelstam variables in the decay process of Y(mS)(pa) —> 
(nS) (pp) P (pe) P(pa) 


s—(pe- pa),  tp=(Pa— Pp), Up=(Pa—Pa)’, 


3sop = s + tp + up = mms) T mA Gus) T 2m2, (7) 
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where P denotes the pseudoscalar m or K, since we also need to take into account the virtual process 


T(mS)(pa) > Y(nS)(py) K (p.)K(p4) in the coupled-channel FSI. £p and up can be expressed in 


terms of s and the helicity angle 0 according to 


1 1 
tp = 5 [Bsop — s + &r(s) cos 6] , up = 5 [3sop — s — &r(s) cost] , 


4m2 
&p(s) = opal? (mis), (ns) s), op=\j\/1— 2 ; (8) 


where @ is defined as the angle between the initial T(mS) and the positive pseudoscalar in the 
rest frame of the PP system, and A(a,b, c) = a? + b? + c? — 2(ab + ac + bc). We define q as the 
3-momentum of the final bottomonium in the rest frame of the initial state with 
lg = 1 o m 2 
q| 2mr(ms) À (MI msp (s): s). (9) 
The calculation of the tree amplitudes is very similar to our previous study of Y (35) decays [8], 
so here we just quote the partial-wave-projected results. Parity and C-parity conservation require 
the pion pair to have even relative angular momentum l. We will only consider the S-wave and 
D-wave components in this study, neglecting the effects of higher partial waves. For the S-wave, 


the amplitudes of the chiral contact term and the Z-exchange term read 


JE esce s aee my Cl (s — 2mp) +2)s+@(1-%)]} (10) 
0 F2 (mS) (nS) P 2 3 , 
~ Zbor 2 /MI mS) Y (nS 
Mig (s) = -S Y mzy Cnnif (8 + lal?) Qol) — lalo? [B3 Qoluri) — vri] } 


Fanna) i=1,2 i 
(11) 


where Cmn, = Cz, t(msy&CZut(nS)s: Yri = (3800 — $ — 2m, )/«(s), and Qo(y) is a Legendre 


function of the second kind, 


LJ dz 1, y+l1 
Qo(y) — sf n = = (12) 


(P;(z) refers to the standard Legendre polynomials). Note again that we consider the Zy-exchange 
diagrams only for the process involving pions. For every heavy particle, namely the bottomonia 
and the Zp states here, a nonrelativistic normalization factor VM has been multiplied to the 
expressions, with M the corresponding mass. The widths of the Z states are neglected in the 
present calculation, since their nominal values are of the order of 10 MeV and thus much smaller 
than the gap between their masses and the T (IS)s threshold. 

For the D-wave, in which 77 scattering is elastic to very good approximation in the energy range 


considered, we only consider the single-channel FSI, and therefore we just give the amplitudes of 
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Type 1(c1)|[B, B, B, B*], |B, B*, B, B*|, [B*, B, B*, B], |B, B*, B*, B*], 
B" B, B*, B*], [B*, B*, B, B*], [B*, B*, B*, B], [B*, B*, B*, B*] 


Type 1(d1)|[B, B, B*, B*], [B, B*, B, B*], [D*, B*, B, B, [B*, B, B*, B, [B, B*, B*, B*], 


H* Hm BB? |, [B*, H*.. Hi, B7]; [B*. B*, B* Hl. [B*, B*, B*, H*]. 


TABLE I: All loops contributing in each diagram class. The mesons are listed as [M1, M2, M3, M4], type 
l(cl) and type 1(d1) refer to the corresponding diagrams in Fig. 1. Two more configurations appear as type 
1(cl) in principle, namely [B, B, B*, B*] and [B*, B, B, B*], however, their contributions to amplitude 444, 


see Eq. (16), vanishes, and hence they are strongly suppressed. Flavor labels are dropped for simplicity. 


the process involving pions, 


x 2 
M" (s) = 3p V P r(msym(ns) colq}, (13) 


^ Zb n 5,/MY (mS) MY (nS 
MPT) = -T o D mas Cms [s lal? — lal"czu£] [Byes - DQolri) ~ Bui) 
TaT i=1,2 
(14) 


Now we briefly discuss the calculation of the box diagrams. There are four intermediate bottom 
mesons in the box diagrams Fig. 1 (c1) and (d1), where we denote the top left one as M1, and the 
others as M2, M3, and M4, in counterclockwise order. The individual contributions are listed in 
Table I, with the pseudoscalar or vector content of [M1, M2, M3, M4] explicitly shown. For the 
Y (mS) — Y (nS)K K processes, some intermediate states can be strange bottom mesons BS, and 
there are four possibilities for each |M1, M2, M3, M4] given above. For simplicity, we do not list 
the combinations of intermediate states in the Y(mS) > Y(nS)K K processes explicitly. 


The general amplitude for the process Y(mS)(pa) + Y (nS)(py)P(p.)P(pq4) reads 


Af (Y(mS) + T(nS)PP) = & ims) A9 (Y (mS) 2 T(nS)PP), (15) 


er (nS) 
and .@)(Y(mS) — Y(nS)PP) can be decomposed as 
A (Y (mS) > Y(nS)PP) = P M4..., (16) 


where we have omitted the remaining terms proportional to tensor structures built from the dif- 
ferent momenta. For the loop amplitude, we have checked that the .4 term is indeed numerically 
dominant, which agrees with the argument that other contractions of the polarization vectors are 


suppressed in the heavy-quark nonrelativistic expansion. So in the following we will only keep the 


0.2 04 06 0.8 1.0 12 02 0.44 0.6 0.8 1.0 12 
s [GeV?] s [GeV?] 


GeV? 


0.10 0.15 0.20 0.25 0.30 0.10 0.15 0.20 0.25 0.30 
s [GeV?] s [GeV?] 


FIG. 3: The S- (left) and D-wave (right) box amplitudes in T(4S) > Y(1S)r*«- (top) and T(49) ^ 


YT(2S)x*«- (bottom). The black solid and red dashed lines denote the real and imaginary parts, respectively. 


terms proportional to €y(ms) ` €r(ns), as we did for the tree amplitude. Details on the analytic 
calculation of the box diagrams are given in Appendix A. 
The partial-wave projection of the loop amplitude for the Y (mS) — Y(nS)PP process can be 


denoted as 
ap ne gJHH(mS) 9JHH(ns) AmpBox{ (s). (17) 


The analytic expressions of AmpBox}’ (s) are very involved, so in Fig. 3 we only plot the numerical 
results for T(4S) > Y(1S,2S)x*x~ in the physical region. Note that the imaginary parts, which 
are due to the on-shell BB intermediate states, are tiny due to the smallness of phase space and 


the fact that the BB pair appears in a relative P-wave. 


D.  Final-state interactions with a dispersive approach, Omnès solution 


There are strong FSIs in the mr system especially in the isoscalar S-wave, which can be taken 


into account model-independently using dispersion theory. Based on the principles of unitarity and 
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analyticity, dispersion theory determines the decay amplitudes up to some subtraction constants, 
which can be fixed by matching to the results of chiral effective theory. Since the mass difference 
between the T(4S) and the Y(1S) is larger than the KK threshold, we will consider the isospin 
symmetric two-channel (mz and KK) FSI for the dominant S-wave component, while for the 
D-wave only single-channel mm FSI will be considered. 

For the T (mS) — Y(nS)r*«- processes, the partial-wave expansion of the amplitude including 


FSI reads 


lage: 


A (s, cos 8) = er(ns) erims) >. [MP (5) + NEF (s) Fi(cosé), (18) 


l 


ll 
o 


where MF (s) represents the right-hand cut part and accounts for s-channel rescattering, and the 
“hat functions" MF (s) contain the left-hand cuts, contributed by crossed-channel pole terms or 
open-flavor loop effects. In general the box diagrams contribute to both the left-hand cuts at t, u > 
(Mp) Mpo)? and right-hand cut at s > (mpo) M pgo) )?, however, this right-hand cut is far away 
from the physical region, so it can be safely neglected. In this study, we approximate the left-hand 
cuts by the sum of the Z;-exchange diagram and the box diagrams, Mf (s) = M; Zo (s )- M, Joop e g), 
Similar methods to approximate the left-hand-cut structures by including resonance exchange (in 
the case of no loops) have been applied in Refs. [28-31]. 

Next we discuss the Omnes solution to obtain the amplitude including FSI. For simplicity first 
we discuss the single-channel solution, which applies for the D-wave case. The functions Mis) are 
real along the right-hand cut, so in the elastic mm rescattering region the partial-wave unitarity 


conditions reads 
Im Mj(s) = [Mi(s) + Mi(s) | sin 89 (s)e- 191 (9) , (19) 


In the elastic region, the phases Of of the partial-wave amplitudes of isospin J and angular momen- 
tum / equal the 77 elastic phase shifts modulo nz, as required by Watson’s theorem [32, 33]. The 
Omnés function is defined as [34] 


f(s) = exp {=f gio, (20) 


T Jamz BOBS 


which obeys Q} (s + ie) = eiði O (s — ie). Using the Omnès function, the solution of Eq. (19) can 
be obtained [8, 35] 


s” [> dz x) sin ôl (x 
mil) = 0) Pr + = ne oars (21) 


T Jamz X" |OF(x)|(z — s) 


10 
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where the polynomial P^ !(s) is a subtraction function. For the D-wave phase shift ó9(s), we will 
use the result given by the Madrid-Kraków collaboration [36], and smoothly continue it to 7 for 
8 — oo. 

For the S-wave, we will take into account the two-channel rescattering effects. Along the right- 


hand cut, the two-channel unitarity conditions reads 
Im Mo(s) = 29" (s)25) [Mo(s) + Mo(s) , (22) 


where the two-dimensional vectors Mo(s) and Mo(s) contain the right-hand- and the left-hand-cut 


parts of both the mz and the KK final states, respectively, 


Mà (s " ME (s 
mo(s)=( 8) | E = ( EO) |, (23) 
FaMo (s) FaMo (s) 
The two-dimensional matrices T$ (s) and X(s) are 
9(¢) 62459 (5) _ . 
T9(s) — ur EE = EOE 4 
0 (s) = O(s)e (4060-30) 1 (2 ) 


a(o) Jeee ai 


2i (s) 
and X(s) = diag(o,(s)6(s — 4m2),o0K«(s)0(s — 4m2.)). There are three input functions in the 
TO(s) matrix: the mm S-wave isoscalar phase shift 69(s), for which we will use the result from the 
Roy equation analysis in Ref. [37]; the 7m > KK S-wave amplitude gg(s) = |gf (s)|e? 99) with 
modulus and phase, for which the results based on the Roy-Steiner approach in Ref. [38] will be 
used. These inputs are used below the appearance of additional inelasticities from 47 intermediate 
states, namely up to \/s9 = 1.3 GeV (the fo(1370) resonance is known to have a significant coupling 
to 4v [39]). Above so, the phases 69 (s) and và are guided smoothly to 2m [40] 
2 
1+ (22° 


$0 


ó(s) = 2m + (ó(so) — 27) 
The inelasticity (s) in Eq. (24) is related to the modulus |g9(s)| by 
mp (s) = y1— 4es(s)ox (s)lgp (s) P0(s — 4m3.). (26) 
'The numerical solution of the homogeneous two-channel unitarity relation 
Im (s) = WT (s)X(s)Q (s), Q(0) 2 1, (27) 
has been computed in Refs. [40-43]. Using Q(s), the solution of the inhomogeneous two-channel 
unitarity condition in Eq. (22) is given by 


s^ [9 dx Q-l(x)T(z)X(x)Mo(x) 
f — (28) 


TL 
m2 T =s 
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To determine the required number of subtractions that guarantees the dispersive integrals in 
Eqs. (21) and (28) to converge, we need to investigate the high-energy behavior of the integrand. 
First it is known that for the single-channel Omnés function defined in Eq. (20), it falls off asymp- 
totically as s^ if the phase shift ôl (s) approaches kr at high energies. Since the D-wave 77 phase 
shift, ó9(s), reaches 7 for high energies, we have Q9(s) ~ 1/s for large s. Second, the high-energy 
behavior of the two-channel Omnès function has been analyzed in Ref. [40], and the 1/s asymptotic 
behavior of Q} (s) is ensured by the asymptotic condition 370; (s) > 27 for s — oo, where » ôf (s) 
is the sum of the eigen phase shifts. 'Third, we have checked that in an intermediate energy range 
of 1 GeV? < s &« mî, both the inhomogeneities contributed by the Z)-exchange term and the box 
graphs term grow at most linearly in s. So we conclude that in the dispersive representations for 
Mo(s) and Mo(s), three subtractions are sufficient to make the dispersive integrals convergent. 

At low energies, the amplitudes M»5(s) and Mo(s) should match to the results of chiral per- 
turbation theory. Namely, in the limit of switching off the final-state interactions, O9(s) = 1 
and (0) = 1, the subtraction functions agree with the chiral representations given in Eqs. (10) 
and (13). Since both MxX(s) and M3(s) grow no faster than ~ s?, they can be covered by the 


degree of the subtractions. Therefore, for the D-wave, the integral equation takes the form 


s? f° dz MZ (zx) sin ó9(x 
uro) = ati) bane +S [^ Soe (29) 


For the S-wave, the integral equation reads 


Mo(s) = Q() (Mj) += f 


T 


* de Q7 (z)T(r)E(z)Mo(z) 


3 E 
m2 T us 


(30) 


where MX(s) = (MX’"(s), 2//3 MX*(s))*. 
The differential decay width for Y(mS) > Y(nS)n*n«- with respect to the v7 invariant mass 


and the helicity angle reads 


dr vsas|q| | ^ " 2 
—— > m Me ME MZ Mr P. 6 3] 
ds d cos 0 1281?m3 (ms) o + Mj + (MZ + Mz) P2(cos 8) (31) 


III. PHENOMENOLOGICAL DISCUSSION 


The experimental data considered in this work are the mr invariant mass distributions of the 
Y (4S) — Y(1S,2S)r+x~ decays measured by the BaBar [44] and Belle Collaborations [45]. 
The chiral coupling constants c; in Eq. (1) are different for the two decays T(4S) > T(1S)r^«z- 


and Y(49) > Y(2S)z*7-, since there is no symmetry connecting the bottomonium states with 
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FIG. 4: Fit results of the 77 invariant mass spectra for the decays Y(4S) > Y(1S)m*«- (left) and T(4S9) > 
Y(2S)n*«- (right). The solid squares and solid circles denote the data from the BaBar Collaboration [44] 
and Belle Collaboration [45], respectively. Fit I (green dash-dot-dotted): only including the chiral contact 
terms c;, Fit II (red dashed): chiral contact terms and the Z,-exchange term, Fit III (blue dot-dashed): 
chiral contact terms and box diagrams, Fit IV (black solid): including the contact terms c;, the Zy-exchange 


term, and the box diagrams. FSI is included in all fits. 


different radial excitations. The mass difference between the two Zp states is much smaller than 
the difference between their masses and the Y(LS)« (l = 1,2) thresholds as well as my(4g) — ma; 
they have the same quantum numbers and thus the same coupling structure as given by Eq. (3). 
So the Z,(10610) and Z,(10650) contributions are strongly correlated in the fit, and it is very 
difficult to distinguish their effects from each other in the processes studied in this work. Therefore 
we only use one Zp, the Z,(10610), in our fit by setting Cnm,2 = 0 as in our previous analysis 
of Y(3S) + Y(1,2S)«« [8]. For the input mass of the Z,(10610), we will use the heavy-quark 
spin symmetry conserving result given in Ref. [46]. The value of gy 49) is extracted from the 
measured open-bottom decay widths of the T (4S), g;g gr(4s) = 1.43 Gey. 

For each T(4S) + Y(nS)r« (n = 1,2) process, the unknown parameters are c; and co corre- 
sponding to the chiral contact TY’@® coupling, C4,,1 related to the Z, exchange mechanism, and 
9JHH(ns) for the box diagrams. To illustrate the effects of the Zy-exchange and the box graph 
mechanisms, we perform several fits by choosing different strategies. Fit I only includes the chiral 
contact c; terms; Fit II adds the Z-exchange terms to them. Fit III includes the chiral contact c; 
terms and the box diagrams, and finally Fit IV takes all of the contact c; terms, the Z exchange, 
and the box diagrams into account. FSI is included in all fits. In Fig. 4, the fitted results of 
Fits I-IV are shown as the green dash-dot-dotted, red dashed, blue dot-dashed, and black solid 
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TABLE II: Fit parameters from the best fits of the T (4S9) > Y (nS)»« (n = 1,2) processes. 


Y(4S) B Y(19)n*n- Y(4S)o Y(2S9)n*n- 


c [GeV] (9.8 + 1.0) x 1074 (1.2 + 0.6) x 107! 
c2 [GeV !] (-1.66311)x10-* (-—1.0+0.6)x 10-! 
Cani (2.6 ge 1.3) x 1074 (—3.2 E: 1.8) x 1072 


(8.6+6.1) x 1075 (1.7 x 0.8) x 107? 


—3 
gJHH(ns) [GeV ?] 


x?/d.o.f 10.45/(20 — 4) 20.65 — 0.04/(7 — 4) = 0.01 


lines, respectively. The fitted parameters as well as the x? /d.o.f. of our best fit, Fit IV, are shown in 
Table II. We find very different values for the parameters c; and c» from fitting the data of transi- 
tions between different Y (IS) states. These low-energy constants parameterize the nonperturbative 
QCD matrix elements of gluonic operators between the initial and final bottomonia. For different 
initial and final Y states, these parameters are not related to each other at the hadronic level, and 
can well be very different. In principle, the parameter values from the fit in this paper cannot be 
directly compared with those in Ref. [8], which do not include the box diagrams when analyzing 
the Y(3S) and Y(25) dipion transitions. We thus made a new fit to the decay T(3S) > Y(1S)a7 
studied therein. It turns out that the values of cy and cy decrease only by around 35% in compari- 
son with those given in Table I of Ref. [8]. Our fittings turn out to indicate the following hierarchy: 
Ict2!] < |letz?^| S Ietz2!| « |etz?^|. This may be understood from the node structure of the T 
wave functions: for the processes with the same initial Y state, the larger the difference between 
the principal quantum numbers, the smaller the gluonic matrix elements and thus the magnitude 
of the parameters. Note that the total x? value for the transition T(4S) > T(2S)n*«- is very 
low, x?/d.o.f. = 0.01. This small number reflects the observation that the fluctuation in the data 
appears to be significantly smaller than what the error bars allow for, which indicates that they 
might well be dominated by systematics. 

Using the central values of the parameters in the best fit, in Fig. 5 we plot the moduli of the 
S- and D-wave amplitudes from the c; terms, the Z,(10610) state, and the box graphs for the 
processes T(4S) > T(1S)r*«- and T(4S) > T(2S)r*7' , respectively. In addition, in Fig. 6 we 
show the resulting theoretical predictions for the angular distributions. 

As shown in Fig. 4, including the Z,-exchange and the box graph contributions improves the fit 


quality for Y(49) + Y(1S)x*-- only marginally, mainly in the region around 1 GeV. However, 
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FIG. 5: Moduli of the S- (left) and D-wave (right) amplitudes in T(4S) + T(1S)n^-- (top) and T(4S) ^ 
Y(2S)n*«- (bottom). The black solid lines represent our best fit results, while the red dot-dashed, blue 
dashed, and green dotted lines correspond to the contributions from the c; terms, the Z5(10610), and the 


box diagrams, respectively. 


for T(4S) + T(2S)r*^-, the fit quality increases significantly when considering either of those 
two mechanisms (or both). Loop effects were already studied in the 3P) quark-pair-creation model 
in Ref. [47], and found to be tiny for Y(39,2S) > T(2S,1S)n^--. This is probably due to the 
fact that 'T (39,2) are too far below the BB threshold. This situation is expected to change for 
the Y(4S), with the open-bottom channels contributing significantly to its decay rate. In Fig. 5, 
one observes that for the dominant S-wave amplitudes, the contributions from the c; terms, from 
the Zy-exchange term, and from the box diagram term are all of the same order. Especially, for the 
decay Y(A4S) > T(1S)z ^7, the box graphs and the Z exchange play a major role in the energy 
range around 0.95 GeV, and account for the better description of the data there. Note that the 
contribution of loops including B, mesons, producing kaons that subsequently rescatter into a pion 
pair, is entirely negligible: in the NREFT formalism, these graphs vanish at the K K threshold. 


For the D-wave, the contributions from Z, exchange and the box graphs are much smaller than 
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cos 0 


FIG. 6: Theoretical predictions of the helicity angular distributions for the decays T(4S) > Y(1S)r*«- 
(left) and Y(4S) + T(2S)r*7- (right). The line style is as in Fig. 4. 


that from the c; terms. We should mention that the plots in Fig. 5 correspond to using the central 
values of the best fit parameters. The shapes of the curves corresponding to the box diagrams and 
the Z,-exchange terms are similar; however, their relative strength is not very meaningful because 
there is a strong correlation in the fit between the parameters C411 and J;gg(1s). This can be 
easily seen from the fact that the curves for Fit II and Fit III are very similar to each other in 
Fig. 4 (left), which means that the Zy-exchange and box terms can hardly be distinguished in the 
c7 invariant mass distribution of the transition Y(4S) > Y(1S)rta-. 

Notice that in Refs. [48, 49], the loop contribution of the sequential process T(49) > BB — 
Y(nS)S — Y(nS)r*z-(n = 1,2), where the scalar S can correspond to the fo(500) and the 
fo(980), has been considered. This kind of loop topology can be described by Fig. 2 (a) including 
FSIs, which is suppressed compared to the box graphs in NREFT. In our scheme, the FSIs are taken 
into account in a model-independent way, and we do not have to specify the contributing scalar 
resonances. Another merit of our calculation is that, instead of only obtaining the absorptive part 
of the loops by using Cutkosky rules |48, 49], we completely compute both their real and imaginary 
parts. 

An interesting feature of the mr invariant mass distribution of T(4S) > Y(1S)r*-- is that the 
older Belle data from Ref. [50] hint at a two-peak structure in the range of m,+,- = 0.8... 1.2 GeV, 
while the later measurements given in Refs. [44, 45] do not display such a feature in any obvious way. 
As the mass difference between Y(4S) and Y(1S) is about 1.12 GeV, the isoscalar-scalar f9(980) 
meson, which couples strongly to 77, should be visible in the spectrum. With FSI described reliably 
in the dispersive approach, we see that the f9(980) indeed accounts for a dip at its mass, and a 


two-peak structure is naturally produced. A possible reason why such a two-peak structure is not 
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observed in Refs. [44, 45] may be the wide energy bins used in these experimental measurements. 
The fact that the f9(980) should be manifest in the mr invariant mass distribution of T(4S) > 
Y(1S)r*-- has already been emphasized in Ref. [7]. The dip caused by the fo(980) is also present 
in the calculation of Ref. [51]. 

For the Y(4S) > T(2S)r*«- process, it is known that the two-hump behavior in the mz 
invariant mass spectra is incompatible with the prediction from the QCD multipole expansion, 
resembling the case of Y(3S) > Y(1S)z [5, 8, 47]. In the formalism outlined above, the original 
formulation of the QCD multipole expansion appears by including only the tree-level c;-terms, 
however, omitting the 77 FSIs. As shown by the blue dot-dashed line in the right panel of Fig. 4, 
including the final-state interaction can roughly reproduce a two-hump structure. However, it 
produces a zero in the amplitude inside the physical region and the agreement with the data is not 
very convincing. This feature was also observed in our previous study of T(3S) — Y (1S), where, 
however, a simultaneous fit of the 77 invariant mass and the helicity angular distributions cannot 
reproduce the two-hump behavior in the dipion mass spectra by only using the c; terms [8]. The 
angular distribution data are therefore important to distinguish the effects of different mechanisms. 
In Fig. 6, the theoretical predictions of the helicity angular distributions in different fit scenarios are 
shown. For T(4S) > T(2S)r*7 ^, the angular distributions are distinctly different when including 
the Z,-exchange and box graphs terms, hence these results can be used to check their effects when 
experimental data become available in the future. 

Using the fit parameters given in Table IL, we can predict the decay width of T(4S) — 
Y(1S)K* K-, as well as the corresponding KK invariant mass distribution. The relevant Feyn- 
man diagrams can be obtained by replacing all external pions by kaons in Fig. 1, but without 
diagram (b1) due to the absence of a Z,YK vertex. The Z, contributes also to KK through dia- 
gram (b2) due to the final-state interactions that, especially around the K K threshold, provides 
strong mm — KK transitions. Most ingredients of the amplitude of the Y(4S) —^ Y(19)K*^K- 
process have been given in Sec. II. We omit the KK D-wave, which is negligible due to its 
strong near-threshold suppression. Within lo uncertainties, the prediction of the decay width 


of T(4S) > T(19)K* K- is 
Ixasyoraüs)k*k- = 0.18350 keV , (32) 


corresponding to a branching fraction of TEE x 1075, and the dikaon invariant mass spectrum 
is given in Fig. 7 (top left). The rapid rise of the KK invariant mass distribution in the near- 


threshold region is a result of the f9(980), in line with the dip around 1 GeV in Fig. 4. Like the 
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(4S) > Y(1S)x*«- process, there is a strong correlation between the Z)-exchange terms and the 
box diagrams in the T(4S) + Y(1S)K+tK~ process, and in Fig. 7 we also plot the contributions 
from the the c; terms (top right), the Z;(10610) state plus box graphs (bottom left), and their 
interference (bottom right), respectively. One finds that for the central values of the theoretical 
predictions, the Zy-exchange term and the box graphs nearly cancel each other, and the total line 
shape is quite similar to the c; terms only. Both the rapid rise in the myx distribution and the 
nontrivial structure in the large Myr region of the dipion invariant mass distribution are due to 
the final-state interactions between the light mesons, depicted in Fig. 1 (c1, d1,a2,...,d2), which 
receive contributions from both the Z,-exchange and box diagrams. As a result, their strong 
correlation in the fit to the data of the dipion transitions leads to the significant cancellation in 
the prediction of the myx distribution. The large spread mainly comes from the uncertainties 
of Z,(10610) plus box graphs, and the interference term. These predictions encourage future 


experimental measurements in this channel. 


IV. CONCLUSIONS 


We have studied the effects of Zp exchange and bottom meson loops in the decays T(4S) — 
Y(nS)x« (n = 1,2). The bottom meson loops are treated in the NREFT scheme, in which the 
power counting rules indicate that the box diagrams are dominant. The strong FSIs, especially 
the coupled-channel FSI in the S-wave, are taken into account model-independently by using 
dispersion theory. The forms of the subtraction functions are obtained by matching to the leading 
chiral contact terms. Through fitting the data of the mr invariant mass spectra, the couplings of the 
Y Y'zz and YB) B®) vertices, as well as the product of couplings of the ZY r and ZyY'z vertices 
are determined (where Y and Y’ denote the final- and initial-state bottomonia). For the dominant 
S-wave component, it is found that the contributions from Zr exchange, the loops, and the chiral 
contact term are of the same order. For Y(4S) > Y(2S)r*7-, including the Z,-exchange term 
and the bottom meson loops naturally describes the two-hump behavior in the mr invariant mass 
distribution. Unfortunately, the present data are insufficient to distinguish between the effects of 
the Z, exchange and the bottom meson loops. We provide theoretical predictions of the helicity 
angular distributions, which may be useful to identify the effects of Zy-exchange and bottom meson 
loops with future experimental data. For the T(4S) — Y(1S)r*-- decay, we expect that there 
is a dip in the mr spectrum around 1 GeV, caused by the opening of the KK channel near the 


fo(980) resonance. This dip has probably not been observed yet in the present experimental data 
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FIG. 7: Theoretical prediction of the KK invariant mass spectrum for the decay Y(4S9) > T(19)K*K- 
(top left). The contributions from the c; terms (top right), the Z,(10610) state plus box graphs (bottom 
left), and their interference (bottom right) are also depicted. The shaded areas corresponds to the error 


band. 


yet due to lack of sufficiently precise energy resolution. Improved data to resolve this issue is 
eagerly awaited. We also predict the decay width and the KK invariant mass distribution of the 
(4S) > YT(1S) K* K- process, demonstrating the usefulness of this additional measurement that 
should be feasible at Belle II. 
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Appendix A: Remarks on the box diagrams and four-point integrals 


In this appendix, we will discuss the calculation of the amplitudes that involve four-point loop 
integrals in some detail. We will start by discussing the parametrization and simplification of scalar 
four-point integrals. Then we will introduce a tensor reduction scheme to deal with higher-rank 
integrals. Finally, we will give the leading part of the corresponding integrals (proportional to 


ey: er) for the possible intermediate bottom mesons. 


1. Scalar four-point integrals 


Because of the simpler structure we begin with the first topology as shown in Fig. 8. The 
corresponding scalar integral, evaluated for the initial bottomonium at rest (p = (M,0)) and 


labelled J(09 to be consistent with Fig. 1, reads 


Joe) 


if d*l 1 
(21)* [P — mz  ie]((p — 1)? — m3 + ieJ[(E — q1 — 22? — má + ie — q1)? — mi + ie] 


Q 


ed / dl 1 
l6mimamama J (21) L ~ m+ ie i - M+ fo m — ie 
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1 
——— ———— ay (A1) 
liq? : SET. . 
| -qi-qi- cra — m3 + ie |i -qi- nr — ma + ie 


Performing the contour integration is straightforward since only one pole is located in the upper 


half-plane. We find 


| H12H23124 f dà E (A2) 
2mj1m»omama (21)? (1? + C42 — ie] [I2 + 253] - q + C93 — ie] (1? = 2l - q1 + C24 — ie] i 


where we defined 


2 
q 
C12 = 2019 (m4 + m3 — M), on = 3s (mem - M +a die). 
2 
MiM; 
n = 2a (me ma - M e dd e BE), hij — —. (A3) 
2m4 Mi + Mj 


For the second topology we immediately find 


(Od) = į SHOE 0000 
HSS | are Se aU) mirid ai)? m + te 


—i f al 1 
16mymgmgm4 J (21) E ——— — ie G dc m- ie 


—————— ————— (A4) 
| +q-M+ ED + m3 — ie e -qp- Ga -m + ie 


ma 


Here the possibility for two different cuts to go on-shell leads to a slightly more complicated three- 


dimensional integral 


M1234 / dil 1 
( 


2m1m»mama 22)? [P + dig — ie] (I? = 251 -qi = 254] -q2 + d34 — ie] 
H24 H13 
reo | ( ) 
where we defined 
di, % 
di» = 2m2 (m1 + m3 — M), dui = pi (ms tma d o EL). 
M4 2ma 


2 2 
day = 2m (m emi - M+ qf + SE), dis = ina (m ma Medie). (A6) 


In both cases the remaining three-dimensional momentum integration needs to be carried out 


numerically. 


2. Tensor reduction 


Since each of the interactions of an Y with a pair of bottom mesons scales with the momentum 


of the latter we will have to deal with 


— 121231124 / dl fO (A7) 
2mymamsma J (27)? [P + cra — ieJ? + 22231. q + cz — ie] [I — 2241. qn + c24 — ie] 
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Intermediate mesons | Amplitude 


P, P, P, V 8q1 - qo J O28 
PV —8g, + g; J D^ 

V, P, V, P 45 qi (2j vi +q- qJ) 

P, V, V, V 45% (qiq -q2 — dg  q1) JO! 

V, P, V, V 46% qi (2qh v +q- qJ) 

VVP: V 45% (q-q — dáq q) IVE 

Aara 88i qi . nT + 460P gi (2d ja +q: q,JO*) 

EUNT 45% (A5q, -JOY — aa IY! + ad (a aI! — 4g] JO) ) 


TABLE III: All loops contributing to topology (c1) in Fig. 1. The mesons are listed as [M1, M2, M3, M4], P 
and V denote intermediate pseudoscalar and vector mesons, respectively. The different flavors are dropped 


for simplicity—the full amplitude contains the sum of all possible ones. 


where f (I) = (1, I’, U} for the fundamental scalar, vector, and tensor integrals, respectively. Us- 


ing the momentum of the final state Y, q, and qı = q1—4q(q:q1)/q?, a convenient parametrization 


reads 
JOi = Hates Hoa J dko o oo EO oo 
2mımzmama J Qm [P +e — iE — 2731 a + cz — ie] — 2531 qi + es — ie] 
= gis?) i qi Jo? (i 
and 
gid — ZMH 7 dio 0000000000000 BB 0000000000008. 
2mimamama J (27)? [P + e — tel[l? — 253l: a + ea — iell? — 2541 qi + es — ie] 
4 qq dd ae cds ug. oed. uidi dedi d! un: 
= [ 54 eae EE JE? 4 TH yo , AA geo | TE T.U go. (A9) 
M d qi alla. | 


where the scalar integrals a can easily be disentangled and have to be evaluated numerically. 
'The corresponding expressions for topology II can be obtained by changing the denominators 


accordingly. 
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Intermediate mesons| Amplitude 


P,P, V, V 8qı + qo J^ 

P, V, P, V 48 (a ' qo (2633 JC + GJ) -É (2g + a? JO) — 85% 95 qi -qg JOY 
V,V,P,P 46% gi (20f.108 Fi ja) 

MEME 4o*^gi (u -q2 + 29 q2) FE gj 0) 

P: VVV —4ó^ (aid Oe tdi (qı - a2 — q?) Jo) 

V,P.V,V -Ai"^qi (2a). 93 — (q- qa + 2i - a2) J*) 

WVV,PV —4ó^ (2aigh 74 di (m^q — m?) Jo») 

VV VP = 45% ¢) (2qh JG + (g? — qi - q2) JO?) 

V, V, V, V 482b (469a, qo JAI — hq- q JY! + gi (4 . godt — 4g] J) ) 


TABLE IV: All loops contributing to topology (d1) in Fig. 1; see Table III for further notation. 


3. Amplitudes 


Tables III and IV list the relevant amplitudes for this calculation. We will only give the dominant 
amplitudes, i.e. the ones that contribute to the part proportional to e(Y^) - «(Y) as was explained 
in the main text. We further notice that all box diagrams are proportional to the overall factor 
e" (T^) &() g2 oy mu gung F2. 

Finally, we need to consider the different flavors of the intermediate bottom mesons. For 
topology (cl) with a pair of charged pions four possibilities exist: [B()*, B()-, B(9*, pC90]. 
[B(9-, p(9*, p(9-, BH], [BHO BHO BH? B+), and (BO, B®, B®), B)-]. For topology 
(d1) this reduces to just two: [B®+, B()-, B®), BO) and [B@, BG), B+, B()-1.. For the 
case of neutral pions the number of possible diagrams doubles—a factor 2 that is balanced by the 


factor V2 in the SU (3) light-meson matrix. 
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